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' We study by numerical simulations the transition from the metastable "false vac- 

i uum" to the broken symmetry phase in the preheating stage after cosmic inflation 

V^*) , in a scalar hybrid model. We take quantum fluctuations and their back reaction 

■ into account by applying a one-loop bubble- resummation. 1 
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1. Introduction 



At present there is a wide class of inflationary models, having in common 
O |' the aim of explaining various kinds of cosmological puzzles 2 . One of the 

consecutive questions raised is how inflation has ended. In a scenario called 
, "Hybrid inflation" 3 the inflationary expansion ends by a phase transition 

from a metastable "false vacuum" to the true, i.e. broken symmetry vac- 
uum. The subsequent phase of (p)reheating 4 after inflation is characterized 
by a phase of abundant particle production due to spinodal amplification 
5^ ■ of quantum modes and due to parametric resonance. We will study some 

aspects of this phase in the following. 
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2. Effective action and renormalized equations of motion 

The Lagrangian for the Hybrid model is given by 

C = \d^d^ + \d»Xd»X - irn 2 $ 2 - \g 2 $> 2 X 2 - ^(X 2 - v 2 ) 2 . (i) 

The field X has an effective mass square m 2 x ($) = —Aw 2 + <? 2 $ 2 and be- 
comes unstable if the absolute value of <f> is lower than $ c = (spinodal 
region). Both quantum fields $ and X are assumed to have classical expec- 
tation values, i.e. {<£} = <fi (inflaton field) and (X) — \ (symmetry breaking 
field). 

We will use here the so-called two-particle point-irreducible (2PPI) ef- 
fective action (EA) formalism. 5 ' 6 The renormalized equations of motion 
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(EOM) that describe the nonequilibrium dynamics of classical mean fields 
and their quantum fluctuations, all coupled to each other, are derived from 
the EA functional 

T[<t>, x , A W<MX)XX} ] = S[<l>, x\ + y J d D x (A 00 A XX + 2A| x ) (2) 

+T/^ A xx + r2PPI ^x,^ x , xx} ], 
r 2PPI [0,x,^Wxx}] =0 + 0- + ©+M>... ■ 0) 

The classical EOM follow from |^ = and |E = 0, while the gap equations 
are derived from = with ij — (fxf), \X or 4>X- The propagator G fulfills 
a constraint equation 

(G- 1 ) ij (x,x')=i(n5 ij + M 2 ij (x))d^(x-x') . (4) 

It is dressed by a resummation of local self-energy insertions. Assuming 
spatial homogeneous fields, Gy (i, if ; p) can be decomposed in momentum 
space into mode functions via Gy(i,t';p) = J2 a =i 2u~ fj " 1 p) ■ 

The quantities uj a (a = 1,2) have to be determined self consistently 
at the initial time [t = 0). For practical calculations one has to trun- 
cate the infinite series in Eq.(3). Here we will only keep the one- loop 
term which leads to a summation of bubble diagrams. a It is given by 
T^iMl^M^M 2 ^] = §Trln[G- x ]. The tadpole insertions derived from 
denote 

jfp(l) 2 /■ J-D— 1 1 

AS>(*) = = E / (^2^ Re ttr fe*wrfep>] • (5) 

These quantities have quadratic and logarithmic divergencies and have to 
be renormalized. It turns out 1 that the EA is renormalized by a simple 
vacuum counter term 8E wa _ c = -S^M^) 2 + 2{M\ X ) 2 + {M 2 XX ) 2 ], where 
in dimensional regularization 6£ = ~ g^a (f — 7 + 1 + ln47r). The vacuum 
counter term 6E vac indeed reflects a resummation of the standard coun- 
terterms derived from a counterterm Lagrangian. 5 The renormalized gap 
equations for the variational parameters M-lj take the form 

Ml^(t)=m 2 +g 2 [ X 2 (t) + A xx (t)}~Ag 2 6^Ml xx (t) , (6) 
Ml xx (t) = -Xv 2 + g 2 [<j> 2 (t) + A 00 (t)] + 3A [ X 2 (t) + A xx (t)] 

-V^Mi^W - 12\6£Ml xx (t) , (7) 

M 2 KAx (t) = 2g 2 W)x{t) + A 0x (t)] - %g 2 5tM 2 KAx (t) . (8) 



a Equivalcntly one can take the local part of a two-loop-2PI approximation. 
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The final set of gap equations, with all finite parts from dimensional reg- 
ularization, is rather lengthy 1 and we don't write it down here. However, 
one can see from Eqs. (6)- (8) that the gap equations form a (3 x 3) sys- 
tem of linear equations, whose coefficient matrix has to be diagonalized by a 
time-independent rotation matrix. This fixes the various mass and coupling 
constant counterterms, obtained from a standard counterterm Lagrangian, 
in a non-perturbative way. 

After all the classical EOM are given by 

Ht) +Ml <p<p m(t) + M 2 R ^ x (t)x(t) - 2g 2 x 2 (t)<t>(t) = ,(9) 
X(t) +Ml xx (t) X (t) +Ml <px m(t) - 2X X 3 (t) - 2g 2 <p 2 (t) X (t) = 0(10) 

Inserting the decomposition of Gij(t,t;p) into mode functions f"(t,p) in 
Eq. (4) the system of EOM for the fj*(t,p) denotes explicitly as 

f?(t; P ) +P 2 f?(t; P ) + Ml 4>4> (t)^(t;p)+Ml ft>x (t)f^(t;p) = ,(11) 
f«(t-p) + p 2 f x <(t;p) + Ml xx (t)f x *(t;p)+Ml <t>x (t)mt;p) = .(12) 

3. Results 

The EOM (9)-(12) are solved numerically. After rcnormalization the mo- 
mentum integrals are finite and are carried out using a non-equidistant 
momentum discretization. In order to simulate the phase transition at the 
end of inflation we choose the initial amplitudes 0(0) > (f> c and %(0) w 0. 
In particular we take (f>(0) = 1.2 and x(0) = 10~ 7 . The other parameters 
are here m 2 = 0, v = 1, g = 2A, A = 1. 

The time evolution of the classical fields <p and % and of the variational 
masses is displayed in Fig. 1. One observes three time regimes: (i) the 
first regime is characterized by a slow-roll of the field <p, where quantum 
fluctuations are almost negligible; (ii) in the spinodal region the mass M. xx 
becomes negative (Fig. lb). Gaussian fluctuations build up (left column 
in Fig. 2) and drive M xx back to positive values. Along the way the 
amplitude of x growths exponentially; (iii) at later times the momentum 
spectra display several spikes (see e.g. the right column in Fig. 2) and the 
symmetry breaking field oscillates around a nonzero value with a relatively 
constant amplitude. The details in the dynamics depend on the parameters 
in the model while the characteristic features are very similar. 1 

In Fig. 3 we display the time and space evolution of correlations between 
the fluctuations. One can see that spatial correlations build up in the 
spinodal region and propagate by twice the speed of light. The correlations 
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Figure 1. Time evolution of (a) the classical fields (f>(t)/ij> c (solid line) and \{t)/v 
(dashed line), (b) the effective masses Ad 2 j(t) with ij = <f><f> (solid line), ij = XX (dashed 
line) and ij = (f>\ (dotted line); the vertical dotted lines indicate the times 19.5, 51.3 and 



100.2; parameters: g 2 = 2A, <j>(0) = 1.2, x(0) = 1.0 X 10" 7 , 



0, A = 1 and v 1 
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Figure 2. Momentum spectra Gij (t, t; p)p 2 /(27r 2 ) for the simulation in Fig. 1 at the 
times t = 19.5 (left column), t = 51.3 (middle column) and t = 100.2 (right column). 



decrease once the field x starts to oscillate around a non-zero expectation 
value. 




Figure 3. 
Correlation 



function rC xx (r, t) 



r i il^F eip ' r E a 2^ Re [/x/n for thc 

simulation in Fig. 1. The propagation ve- 
locity is r/t = 2 (twice the speed of light) . 



The decoherence time tdec(p) for which a given mode becomes classical 



is defined by \Fij(td ec , p) = 1 with Fij(t,p) = Im£] Q \ u 1 ]■ The modes 
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right to the curves in Fig. 4 never become classical. The functional de- 
pendence is not proportional to p 2 , as found in models with a slow quench 
where M 2 Jt) oc (t - t). 



Figure 4. Decoherence time idec(p) f° r 
which a given mode p becomes "clas- 
sical" for the simulation in Fig. 1; the 
solid line represents the inflaton and the 
dashed line the Higgs modes; the dotted 
line corresponds to t^ cc oc p 1 . 



We have analyzed by numerical simulations the transition from the 
mctastable phase to the broken symmetry phase in the Hybrid model. We 
observed that the back reaction of quantum fluctuations plays an impor- 
tant role and leads, e.g., to a dynamical stabilization in the spinodal region. 
The approximation used here provides reliable information on the nonequi- 
librium evolution at earlier times. However, one has to improve further in 
order to cure, e.g. the lack of dissipation at later times. While there are 
various improvements, once scattering is taken into account via non-local 
two- loop diagrams, even in the 2PPI resummation scheme 6 , one ultimately 
wants to achieve a universal late time behavior, i.e. reaching quantum ther- 
malization. The full treatment of the reheating phase remains a desirable 
goal for future investigations. 
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4. Conclusions 



